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Abstract 

The charge-symmetry-breaking amphtudes for the recently observed dd — > an^ reaction are 
investigated. Chiral perturbation theory is used to classify and identify the leading-order terms. 
Specific forms of the related one- and two-body tree level diagrams are derived. As a first step 
toward a full calculation, a few tree-level two-body diagrams are evaluated at each considered 
order, using a simplified set of d and a wave functions and a plane- wave approximation for the 
initial dd state. The leading-order pion-exchange term is shown to be suppressed in this model 
because of poor overlap of the initial and final states. The higher-order one-body and short- 
range (heavy-meson-exchange) amplitudes provide better matching between the initial and final 
states and therefore contribute significantly and coherently to the cross section. The consequences 
this might have for a full calculation, with realistic wave functions and a more complete set of 
amplitudes, are discussed. 

PACS numbers: ll.30.Hv, 25.10.+S, 25.45.-z 
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I. INTRODUCTION 



For most practical purposes, hadronic isospin states can be considered as charge sym- 
metric, i.e., invariant under a rotation by 180° around the 2-axis in isospin space. Charge 
symmetry CS is thus a subset of the general isospin symmetry, charge independence CI, 
which requires invariance under any rotation in isospin space. In quantum chromodynamics 
QCD, CS means that the dynamics are unchanged under the exchange of the up and down 
quarks Q]. In the language of hadrons, this symmetry translates into, e.g., the invariance 
of the strong interaction under the exchange of protons and neutrons. However, since the 
up and down quarks do have different masses {rriu ^ m^) , the QCD Lagrangian is not 
charge symmetric and neither is the strong interaction of hadrons. This symmetry violation 
is called charge symmetry breaking CSB. There is also a contribution to CSB because of the 
different electromagnetic interactions of the up and down quarks. 

Observing the effects of CSB interactions therefore provides a probe of m„ and m^, 
which are fundamental, but poorly known, parameters of the standard model. The quantity 
rrid is larger than m„, causing a specific pattern of mass splitting between members of an 
isospin multiplet P| . In particular, the light quark mass difference causes the neutron to be 
heavier than the proton. If this were not the case, our universe would be very different, as 
a consequence of the dependence of Big-Bang nucleosynthesis on the relative abundances of 
protons and neutrons. 

Expen^enu, evdence fo. CSB .ee„ den,o....ed . .ix,„, g nucleon 

mass splitting, the binding-energy difference of mirror nuclei such as and ^He Q], the 
different scattering lengths of elastic nn and pp scattering ^, and in the minute but well- 
measured difference between the proton and neutron analyzing powers of elastic np scatter- 
ing 01 . A recent theoretical analysis of vrA^ scattering data found a small CSB effect jsl. 

Studying the dd —>■ an^ reaction presents exciting new opportunities for developing the 
understanding of CSB. This reaction obviously violates isospin conservation; but more specif- 
ically, it violates charge symmetry since the deuterons and the a-particle are self-conjugate 
under the charge-symmetry operator, with a positive eigenvalue, while the neutral pion wave 
function changes sign. This reaction could not occur if charge symmetry were conserved, and 
the cross section is proportional to the square of the CSB amplitude. This is unique because 
all other observations of CSB involve interferences with charge symmetric amplitudes. Thus 
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a very clean signal for CSB is obtained through the observation of a non-zero cross section. 
Furthermore this process has a close connection with QCD because chiral symmetry plays 
a dominant role in determining pion-production cross sections. 

Lapidus, in 1956 9|, was the first to realize that the dd a;7r° reaction would be a useful 
probe of CSB. Various experimental groups tried to observe it, but without success 



After other attempts yielding only upper limits [ll|, a group at the Saturne accelerator in 
Saclay reported a non- vanishing dd — > avr^ cross section at = 1.1 GeV 12] ■ This finding 
was refuted by members of the same collaboration who argued that the putative signal for 7r° 
production actually was caused by the dd 077 background 3]- The importance of this 
background was confirmed by calculations of the double radiative capture 1J| , using a model 
based on a very successful treatment of the dd anir reaction at similar energies [l5| . Thus 
the Saclay dd avr^ cross section is almost certainly a misinterpretation of a heavily-cut 
smooth dd 0:77 background [14 1 . 

There have been two exciting recent observations of CSB in experiments involving the 
production of neutral pions. Many years of effort have led to the observation of CSB in 
np dn^ at TRIUMF. After a careful treatment of systematic errors, the CSB forward- 
backward asymmetry of the differential cross section was found to be Af^ = [17.2±8(stat) ± 
5.5(sys)] X 10~^ [ifjj. In addition, the final experiment at the lUCF Cooler ring has reported 
a very convincing dd —>■ air^ signal near threshold {a = 12.7 ± 2.2 pb at = 228.5 MeV 
and 15.1 ± 3.1 pb at 231.8 MeV), superimposed on a smooth dd 077 background [l7|. 
This background is roughly a factor two larger than calculations based on Ref . , but has 
the expected shape. The data are consistent with the pion being produced in an s-wave, as 
expected from the proximity of the threshold (T^ = 225.6 MeV). 

Clearly, these new high-quality CSB experiments demand a theoretical interpretation 
using fundamental CSB mechanisms. At momenta comparable to the pion mass, Q ~ m^r, 
QCD and its symmetries (and in particular CSB) can be described by a hadronic effective 
field theory EFT, chiral perturbation theory xPT jisl. Ilfll|. This EFT has been extended to 
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pion production [2ij|, |21|, |2^ l23L l24l| where typical momenta are Q ~ v^m^M, with M the 
nucleon mass. (See also Ref. [25[ where pion production was studied neglecting this large 
momentum in power counting.) This formalism provides specific CSB effects in addition to 
the nucleon mass difference. In particular, there are two pion-nucleon seagull interactions 
related by chiral symmetry to the quark-mass and electromagnetic contributions to the 
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nucleoli mass difference l26l hi I . 

n n 

It was demonstrated for the CI reactions tttt tttt M, ttN ttN nj, and A^A^ ^ 
A^A^ [29] that the values of the low energy constants can be understood as the low energy 
limit of the exchange of a heavy state. This procedure is called the resonance saturation 
hypothesis. Within this scheme the other CSB interactions, also caused by the light quark 
mass difference can be viewed as the low-momentum limit of standard meson- 

exchange mechanisms, such as n-rj-r]' and p-u mixing. Determining the various interaction 
strengths may provide significant information about the quark mass difference. Since these 
terms contribute to CSB in the reactions np — > rfvr" and dd an^ with different weights, it 
is important to analyze both processes using the same framework. 



Early calculations of CSB in np dn^ 



32| incorporated most of the relevant mech- 
anisms, giving an asymmetry — dominated by n-r] mixing — of the order of —2 x 10~^ for 
energies near threshold js^. The combined pion-nucleon seagull interactions required by 
chiral symmetry generate a larger contribution with the opposite sign and provide a 
prediction for Af\^{np dn^) (based on a crude estimate of the strength of the CSB rescat- 
tering contribution) that was confirmed by the recent experimental observation. However, 
the experimental value is in the lower band of the predicted range of values of Afb. 

Our aim here is to provide the first study of CSB in the near threshold dd an^ reaction 
using chiral EFT techniques. The effect of -K-rj-r]' mixing on this reaction was studied several 
years ago at = 1.95 GeV 2^. Pion production was assumed to be dominated by the 
production of 77 and t]', followed by tt-t] or Tc-rj' mixing. Using phenomenological information 
on these parameters and on the rj-rj' angle, the cross section was expressed in terms of 
existing data for the 77 production cross sections. This method cannot be used for energies 
lower than that required to produce an t] meson, and other CSB contributions cannot be 
evaluated this way. 

It is necessary to explicitly account for the detailed dynamics of the few-nucleon pion- 
production amplitudes. Therefore we will discuss the CSB amplitudes in the first few orders, 
defined according to a chiral counting scheme that provides a general guide to the expected 
importance of different interaction terms. Such schemes do not explicitly account for spin- 
isospin factors, for the sometimes poor overlap of wave functions, or for the spin and isospin 
dependence of the wave functions. We shall see that selection rules resulting from the use 
of specific wave functions and the threshold kinematics have a strong impact on the relative 
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importance of particular diagrams. 

The fast incoming deuterons (p ~ 460 MeV/c in the center-of-momentum frame cm.) 
need to be slowed down to produce an a-particle and an s-wave pion at threshold. The 
resulting large momentum transfer can be transmitted through the initial- and final-state 
interactions or wave function distortions, and through the exchange of a particle in the pion- 
production sub-amplitude. Only the latter two possibilities will be considered here. The 
complexities of the dd initial state interaction will be included in a future study. Thus, we 
expect that a pion-production sub-amplitude should preferentially provide for momentum 
sharing between the deuterons, in order to avoid forcing the nucleons out into the small, 
high-momentum tail of the a-particlc wave function. 

Spin, isospin, and symmetry requirements restrict the partial waves allowed for the dd — > 
Q;7r° reaction. In the spectroscopic notation ^^~^^Ljl, where S, L, J are the spin, orbital, 
and total angular momenta of the dd state and I is the pion angular momentum, the lowest 
partial waves are ^Pqs and ^Dip. Hence, production of an s-wave pion requires that the 
initial deuterons be in a relative P-wave, with spins coupled to a spin-1 state, coupled 
together to zero total angular momentum. The deuteron spins then need to be flipped, 
while absorbing the P-wave, to form the spin-0 state of the helium nucleus. The invariant 
amplitude therefore takes the form p • (eiXe2) where p is the deuteron relative momentum 
and ei_2 are the polarization vectors of the initial deuterons. On the other hand, a p-wave 
pion is produced only when the deuterons are in a relative D-wave, with spins maximally 
aligned to spin 2, requiring either a couphng with AL = AS — 2 or £)-states of d or a. This 
invariant amplitude is of the form p • eip ■ (e2 x p^) + p ■ e2p ■ (ei x p^), where is the 
pion momentum. Interferences between s and p-waves will disappear for any unpolarized 
observable. 

In addition to these momentum-sharing and overall symmetry considerations, the spin- 
isospin symmetries of the nucleons in the dd : a system will turn out to be crucial in deter- 
mining which sub-amplitudes can contribute and what possible meson exchanges can take 
place. This will be discussed in considerable detail below. 

In this first stage we explore the dd — > avr^ production process using chiral EFT with 
the simplest deuteron and a-particle wave functions, and ignoring the effects of initial-state 
interactions. This will give us an initial test of the amplitudes and provide us with the 
framework necessary to establish the ingredients for a full-fledged model. We are developing 
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a full model, using realistic wave functions and incorporating initial-state interactions, along 
with A admixtures, and the results will be reported in forthcoming papers. 

The chiral power counting scheme is developed in Sec. IHl resulting in a list of possible 
CSB amplitudes. Our simplified model is presented in Sec. Illll The relative importance of 
the amplitudes in this model is investigated in Sec. IIVI The paper then concludes in Sec. |V| 
with a discussion of the results, implications for the interpretation of the lUCF experiment, 
and future prospects. Some details of the calculation are included in an Appendix. 



II. CSB OPERATORS 



We use the EFT power-counting scheme to classify the CSB pion production operators 
in this section. In addition, the specific forms of the tree-level one- and two-body operators 
are derived. A few unknown low-energy constants LECs appear in the first few orders. Since 
these cannot be determined by symmetry considerations, we use phenomenological transition 
amplitudes to estimate their size. The effects of the derived operators are evaluated using a 
simplified model in Sec. IIIII This allows us to check that the leading non- vanishing operators 
of the chiral expansion indeed lead to a CSB cross section of the observed order of magnitude. 



A. Effective Interactions 



In QCD, the pseudo-Goldstone bosons of spontaneously broken chiral symmetry, SU{2) x 
SU{2) — s> SU{2), can be identified with the pions. Chiral symmetry then strongly constrains 
the interactions allowed for pions with matter, and it is possible to construct a well-defined, 
convergent effective field theory for near-threshold pion reactions, namely chiral pertur- 
bation theory. Reviews with special emphasis on nucleon systems are provided in, e.g., 
Refs. The chiral expansion can be adapted to the larger momentum scale inherent 
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23 24|. 
d,|24|. 



in pion production in nucleon-nucleon and nucleus-nucleus collisions 
The necessary power series may converge (albeit slowly) for this class of reactions 
Studies of the pp ppn^ reaction have shown that the resonance-saturation hypothesis does 
not necessarily lead to couplings of natural size, at least for interactions that contribute to 
the production of s-wave pions [21]. This issue should be further investigated. 

We intend to reproduce the S'-matrix elements of QCD at momenta much smaller than 
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the chiral-symmetry-breaking scale, here identified for simphcity with the nucleon mass M. 
To do this, the low-energy EFT must contain all the interactions among pions tt, nucleons 
A^, and Delta-isobars A, that are allowed by the symmetries of QCD. For the following, the 
relevant CI interactions are 



'CI 



1 



iVt[T-(7r X7T)]Ar+|^ 



1 



iVV • a ■ NiV-K) - — iN^T ■na-VN + h.c. 

^ ' 2M 



N^T ■ S ■ A(V7r) + h.c. - — liN^T • tt^ • VA + h.a 



(1) 



Here the first interaction is the Weinberg- Tomozawa term whose strength is fixed by chiral 
symmetry in terms of the pion decay constant /^r = 92.4 MeV. The other terms represent 
the standard axial-vector couplings — including recoil — of the pion to the nucleon (with 
qa — 1.26) and to the Delta-isobar (with Ha — 2.8). Note that a and r are the usual Pauli 

— * 

matrices in spin and isospin space, and S and T are the standard A^A spin and isospin 
transition matrices, normalized such that SiSj' = ^{26ij — ieijkCk) ,TaT^ = ^{26ab — iSabcTc)- 
Charge symmetry breaking can occur either via exchange of a long-wavelength (soft) vir- 
tual photon or via short-range interactions. The former is generated by writing all allowed 
gauge-invariant interactions of the photon field. The latter are represented by local interac- 
tions that come either from the quark mass difference — = e(m„ -|- m^), or from the 
exchange of short- wavelength (hard) photons ("indirect" electromagnetic effects), or both. 
The relevant CSB interactions are 



-^CSB — 



SM 



3SM 

"up 

36M 
'8M2 
1 



T3 - 



TTsT • TT 



T3 - 



TTsT ■ TT 



V^A^+(V^A^)^ (t3- 



N 



N 



TTsT ■ TV — TV Ts 



V^N+{\/^Ny r3 + 



TTsT ■ TV — TV T3 



N 



1 . 



N'^ [-(5MV^ (7r3T • tt) (5M (t^^t ■ tv - 77^3) 



N 



-6M{diNy {'K^T ■ Tv) akdjN + 5M{diN)^ (ti^t ■ tv - 77^3) aud^N 
1 



2M 



iN'^'k^a ■ VN h.c. 



+ ... 



(2) 



where 5M — 0{em'^/M) and 5M — 0{aM/7r) are, respectively, the quark- mass-difference 
and electromagnetic contributions to the nucleon mass difference, and Pi — 0{em1/M'^) and 
= 0{a/7r) are, respectively, the quark- mass-difference and electromagnetic contributions 



to the isos pin - violating pion-nucleon coupling. This Lagrangian is consistent with the one 



from Ref. j2 
from Ref. 13 



with the 5M term added from Ref. [33[ and the pion-nuc 



eon + /^s) term 



An apparent 



This Lagrangian is also consistent with that of Ref. 
difference of an overall minus sign arises because Ref. js^ used different signs for the pion 
field and for 5M + 5M. The CSB seagull term is consistent with the one used in Ref. |33|J . 
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These and other CSB EFT interactions were considered in Refs. 

As usual, we have used naive dimensional analysis to estimate the strengths of the 
terms in the Lagrangian, i.e., we have assumed that the LECs are of natural size. In principle, 
these parameters should be determined using experimental data. We now discuss some of 
the information we have about them. 

The first two terms of Eq. (j2)) are the pion-nucleon seagull interactions required by chiral 
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and can be described as the CSB components of the pion-nucleon a-term. 
The strengths are determined by the coefficients 5M and 5M, with their sum related to the 
nucleon mass splitting: to this order, 

5M + 5M = AM = Mn-Mp = 1.29 MeV. (3) 

The coefficients are not well-known separately. With some assumptions about higher-energy 
physics, the Cottingham sum rule can be used to give 5M = —(0.76 ± 0.30) MeV j35|. 
It is desirable to determine these parameters without these assumptions. The (5M, 5M 
contribution to other observables generally depends on a different combination than that in 
Eq. (jni). It is difficult to isolate the parameters in vrA^ scattering, so it was suggested 
that CSB in pion production could be used instead. The forward-backward asymmetry in 
np dn^ was shown to be sensitive to 6M — 6M/2, but it also depends significantly on 

The other LECs are not well-known either. The pion-nucleon CSB parameter 0-t + 
is constrained by the Niimegen phase-shift analysis of the A^A^ scattering data j36[ to be 
A +/?3 = (0 ± 9) X 10~^ 2D|- Below we estimate the impact of this interaction following the 
standard practice of neglecting and modeling (3i by vr-r^ mixing js^, which is consistent 
with the bound from A^A^ scattering. 

Among the ". . ." in Eq. (j2)) there are several CSB short-range pion-two-nucleon interac- 
tions that contribute in the order we will be considering. One example is 

-^^^^ATtiV lATtaAT.-^ ^ 



2M 



iN^-n^ia ■ VN + h.c. 



(4) 



where we expect that 71 = 0[em^/(/^M^)] and 73 = 0[a/(7r/^M)], for the quark-mass- 
difference and electromagnetic contributions respectively. We know very little about the 
LECs appearing in these short-range pion-two-nucleon interactions, and therefore will model 
these LECs with various heavy-meson-exchange HME mechanisms as detailed below. 



B. Power Counting 

It is necessary to order the various amplitudes according to the size of their expected 
contributions to pion production. There are several strong-interaction scales in the problem, 
namely, 

• X = ~ \Jt^it/M^ the initial cm. momentum of the deuteron divided by the 
chiral-symmetry-breaking scale (here identified with the nucleon mass M), which we 
will use as the expansion parameter; 

• rriT^/M ~ x^i where m^r denotes the pion mass; 



• (j^A ~ M)/M ~ Xi with Ma the Delta mass 5l| — the order assignment given is in 
line with Ref. |[24]; and 

• 7/M ~ x^, where 7 is the typical nucleon momentum inside the deuteron and the a 
particle (for simplicity we will not distinguish between the two). 

Moreover, the strengths of CSB effects are governed by 

• a/vr, the fine structure constant that appears with every exchange of a virtual photon, 
typically with an extra factor of vr; and 

• em^/M^, the factors of — that come from explicit chiral symmetry breaking 
via quark-mass terms [5^ . 

We discuss the two types of contributions individually, to first order, in the following subsec- 
tions. Second-order effects in a and e can also be treated, but are truly small, and ignored 
here. 

Power counting in systems of two or more nucleons is complicated by the fact that some 
diagrams contain small energy denominators, corresponding to states that differ from initial 
and/or final states only by an energy 0(7^/M). Sub-diagrams that do not contain such 
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enhancements are denoted as irreducible. Conservation of energy and momentum in pion 
production requires that at least one interaction takes place among nucleons — before, 
during, or after the pion emission. This interaction transfers a momentum of order p ~ 
y/rriT^M. When such interactions happen before or after pion emission, they are included in 
the (high-momentum tail of the) initial- or final-state wave function. In this case we can 
speak of a "one-body" pion-emission operator. However, in order to compare sub-amplitudes 
of the same dimensions and count powers of we include these interactions as part of the 
irreducible pion sub-amplitude. The full pion-production amphtude is "reducible", because 
it includes further initial- and final-state interactions (via the deuteron and a wave functions) 
that transfer momenta of order 7. 

The separation of reducible and irreducible sub-amplitudes is convenient because it iso- 
lates interactions involving the scale x the irreducible part. Power counting for the 
initial- and final-state interactions corresponding to momenta of 0(7) can be done in the 
usual way In this first paper, we use simple wave functions in lieu of wave functions 
obtained in EFT. The needed EFT wave functions may soon be a reality, since chiral three- 
and four-nucleon calculations already exist 

The loop integrals, propagators and vertices bring factors of momenta, masses, and cou- 
pling constants to any given diagram. Dimensional analysis can be used to express any 
coupling constant as appropriate powers of M times numbers of order 1 (for CI operators) 
or em'l/M'^ 01 a/n (for CSB operators). Some factors, common to all diagrams, are not 
written explicitly. For example, since we study a system of four nucleons that are bound 
in an a particle in the final state, there are always three loops that are controlled by 7. 
Thus, all we need to keep explicitly for a 2n-nucleon operator (in addition to what can be 
read from the vertices and propagators directly) is a common factor (pV(47r)2){"-i) (here 
we have only a three-dimensional integral because we estimate the measure of a convolution 
integral with a wave function). Therefore explicit factors of 7 are not included explicitly in 
the assignments of chiral order. 

A. ...ed .„ Ret. Q Q. .he Me.^cKy of d.a,... . v.. «e.e„t fo. .wave p.o. 
and p-wave pions. We here specialize to s-wave pion production, relevant for the recent 
lUCF experiment. 
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FIG. 1: Leading order diagram with strong CSB. The cross indicates the occurrence of CSB. The 
dot represents a leading-order CI vertex. 



1. Diagrams proportional to e 

At leading order LO there is only one contribution: pion rescattering, where the CSB 
occurs through the seagull pion-nucleon terms linked to the nucleon mass splitting — see 
Fig. m in which the leading CI interaction is represented by a dot, and CSB by a cross. The 
irreducible part of this diagram is 0[em^/(/^Mp)]. The analogous diagram was identified 
in Ref. 3|, using the present counting scheme, as giving the dominant contribution to the 
forward-backward asymmetry in np —>■ dir^. We shall show that, in the dd induced CSB 
reaction, selection rules tend to suppress the rescattering via these seagull terms, if initial 
state interactions are ignored. 

There is no next-to-leading order NLO contribution (suppressed by just one power of 
x). At NNLO, however, there are several contributions, displayed in Fig. |21 The encircled 
vertices stem from sub-leading Lagrange densities. For example, the sub-leading vertex in 
diagrams (a) and (b) arises from the recoil correction of the CSB tiNN vertex, the one in 
diagram (c) denotes the recoil correction of the CI tcNN vertex, and that in diagram (d) 
represents the recoil corrections to the CSB seagulls. Diagram (b) involves the Weinberg- 
Tomozawa vertex. 

Note that diagram (a) can be interpreted as the sandwich of a one-body CSB operator 
between CI initial- and final-state wave functions. It is necessary to include the effects of 
CSB in the wave functions in addition to the diagrams shown in Fig. |21 The easiest way 
to see this is to compare the size of the LO CSB production operator (rescattering via the 
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(a) (b) (c) (d) 




(e) (f) (g) (h) 




(i) a) (k) 

FIG. 2: NNLO diagrams with strong CSB. Vertices with an additional circle originate from sub- 
leading Lagrange densities. We do not display all possible orderings. 

seagull terms) times the LO CI contribution to the A^A^ potential (e.g., one-pion exchange) 
with the LO CI production operator (rescattering via the Weinberg- Tomozawa term) times 
the LO CSB contribution to the A^A^ scattering — assumed to be one-pion exchange with a 
CSB coupling on one vertex. This shows that CSB in the wave functions should be significant 
in a NNLO calculation. Typical diagrams are shown in Fig. El 

The effects of parity conservation suppress the influence of CSB in a single deuteron wave 
function, but CSB does occur in the interactions between the deuterons. One such term 
arises from photon exchange as in Fig. 01 The dominant CSB contribution in the a-particle 
wave function may be expressible in terms of the point radius difference of the neutron and 
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FIG. 3: The influence of strong and electromagnetic CSB in the initial and final state. The 
wiggly line represents the exchange of a photon, while the dashed line represents a meson exchange 
contribution with one CSB vertex. 



proton r„ — Tp, which can be calculated in microscopic models for few-body systems. Results 
of these calculations will be presented in future work. 

Loop diagrams appear already at NNLO. We display only the topology of these diagrams, 
but it is clearly necessary to include all other orderings. A striking feature of the present 
analysis is that, at this order, no counterterms are allowed by the symmetries. The corre- 
sponding counterterms — the CSB four-nucleon contact interactions in Eq. (j2|), displayed 
below in Fig. Hfb), appear first at N^LO. Therefore, those parts of the loops that appear at 
NNLO are to be finite. This situation is in complete analogy to the CI pion production in 
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nucleon-nucleon collisions discussed in detail in Ref. 

Fig. m displays some of the higher-order contributions. A contribution with an intermedi- 
ate A-isobar, that appears at N^LO, is shown in diagram (a). The CSB contact interactions 
displayed in diagram (b) start to contribute at N^LO. Their values will be estimated below 
using phenomenological input. 
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(a) (b) 

FIG. 4: Some typical higher-order diagrams with strong CSB. A double line represents a A-isobar. 
Diagram a) appears at N^LO whereas diagram b) is a N^LO contribution. 

2. Diagrams proportional to a 

Electromagnetic contributions can be ordered relative to each other in exactly the same 
fashion. In this case, the LO is 0[aM/{4:7rf^p)]. These dia grams contain Coulomb interac- 
tions in the initial- or final-state. In particular, the effects of photon exchange between the 
initial deuterons, followed by production by a strong interaction, could be very important. 
An example of such a term is provided by Fig. El 

The NLO electromagnetic diagrams — suppressed by one power of x — that contribute to 
CSB in the production operator are shown in Fig. |S| It is important to note that in threshold 
kinematics (on the two-body level the outgoing nucleons as well as the produced pion are 
at rest) the two diagrams (b) and (c) cancel — in a realistic calculation we should expect 
some of this cancellation effect to survive. The three-body diagram (a) should therefore be 
the one to estimate the photon effects in the production operator at this order. In addition, 
higher-order photon couplings in the wave functions contribute at this order. 

There are various other contributions at NNLO — see Fig. IHl In what follows we will 
explicitly calculate the two-body operator that involves a photon exchange stemming from 
gauging the recoil correction to the ttNN vertex 0|, diagram (a). This will give us an idea 
of the relative importance of soft photons compared to the strong CSB effects. 
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(a) 



(b) 



(c) 



FIG. 5: NLO diagrams with CSB stemming from soft photons. 






(a) 



(b) 



(c) 



FIG. 6: NNLO diagrams with CSB stemming from soft photons. 
C. Heavy-Meson Interactions 

We assume that EFT LECs can be determined using the exchange of massive resonances 
to estimate the impact of short-ran ge p hysics. Such an approach was used in CI pion 



production, for example, in Refs. 
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2l|. In principle the counterterms can be determined 



by other data, and this would eliminate the need for our heavy-meson model. In the present 
context, we include the exchanges of the (a, a;, and p) mesons depicted in Fig. [3 
The meson-exchange diagrams can be calculated from the following Lagrangian: 



HME 



i'- (5) 



Here ip is the Dirac four-component nucleon field and t], a, u^, are the meson fields. We 
use the parameter values in Table H] as representative of typical one-boson exchange OBE 
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n 




0,00,... 




+ 



a) 



b) 



FIG. 7: Resonance saturation for (a) the CSB ttNN vertex modeled here by tt-t] mixing and (b) 
the CSB four-nucleon operators. The elhpsis indicate that additional short-range mechanisms are 
to be included, as discussed in the text. 



models and the standard value Cp = 6.1 for the large ratio of tensor a^j^^d^ /{2M) to 
vector 7^ couphng for the p meson. The //-nucleon coupling g^^ will be discussed below. 



TABLE I: Table of meson masses and coupling constants. 



UJ 



m (MeV/c2) 
550 
783 
770 



SixNN 
Ait 



7.1 

10.6 

0.43 



The photon-nucleon coupling is described by the Lagrangian (up to dimension 5) 



-eip 



1 + T' 



1,A^ + 



2M 



(6) 



where Ao,i = Ap±A„ and Ap = 1.793 and A„ = —1.913 are the proton and neutron anomalous 
magnetic moments. 



D. Explicit form of leading tree- level operators 

We now turn to the explicit form of the leading tree-level two-body operators, in order to 
exploit the selection rules. Corresponding expressions for the loops as well as the three-body 
electromagnetic term mentioned above will be presented in a subsequent publication. 

We start with the formally leading mechanism. Fig. Q together with the recoil correction 
at the pion-nucleon vertex, Fig. |^c). The pion-exchange operator coming from the seagull 



17 



terms is 



r „o 



where = rj — is the relative coordinate of nucleons i and j, kj = — z Vi (k^ = ^ Vi) is 
the initial (final) momentum of nucleon z, = k. — k, is the momentum transfer to nucleon 
i (here symmetrized with the Yukawa factor), and the Yukawa parameter /i = ^J^niT^. In 
our numerical estimates below, we use the value for 6M from the Cottingham sum rule, 
which translates into SM — SM/2 = 2.4 MeV js^. In the fixed kinematic s ap proximation 
for pion production by two nucleons, the exchange pion energy = |39l |. 

It may be noted that the term from Eq. (jHI), proportional to q^, actually gives rise to most 
of the CSB s-wave amplitude in np dii^ js^ . This interferes with CI p-wave production. 
On the other hand, the CSB p-wave amplitude, arising mainly from the CSB one-body 
operator shown in Eq. Q or from a CI production operator following a CSB initial state 
interaction, interferes with the CI s-wave and was about as important in Ref. |33||], but would 
be relatively irrelevant here in the absence of such an important interference at threshold. 

The nucleon recoil term ~ |(k- + kj) is smaller, since it is suppressed by an additional 
factor m-Tr/M. However, if the simple deuteron and a wave functions of Sec. II I II are used, 
the spin-isospin symmetries prohibit this amplitude for nucleons from different deuterons. 
The qj term will integrate to zero inside a single deuteron, leaving the (in-deuteron) recoil 
as the only allowed contribution. Thus the symmetries in this particular model suppress 
the contribution from Fig. ^ leaving only Fig. |21^c): the seagull amplitude is reduced from 
LO to NNLO and there is no momentum sharing. This suppression is expected to be less 
important once initial state interactions are included and realistic wave functions are used. 

At NNLO there are various other contributions. The one-body operator. Fig. |2fa), is 

Oi = ^E'^^-(q^-^(k:+k.))^Aii5:^,.(k:+k.), (9) 

Ai = (10) 

The p-wave q^ = — term is suppressed in the threshold regime considered. In addition, it 
is not allowed in our plane wave approximation, since it lacks the tensor coupling required 
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for the ^Dip transition. The s-wave recoil term is allowed, albeit suppressed by a factor 
uj/M, hence the parameter Ai. This s-wave term is NNLO. 

The isospin- violating (3i is here modeled js^ by tt-t] mixing [see Fig. [7|^a)], 

(3, = -g,{AH\ri)/m% (11) 

where Qr, = grjfn/M = 0.25 is the r]NN coupling constant and {ir^lHlr]) = -4200 MeV^ the 
vr-r^-mixing matrix element 0|- The value of (jr^ corres pon ds to (7^/471 = 0.51, similar to the 
small values found from photo-production experiments |4jJ . However, other values, based on 



hadronic experiments, are as high as g^/4:iT = 3.68 42] or 2-7 for the OBE parameterizations 
of the Bonn potentials Q]. The CD-Bonn OBE potential assumes a vanishing value fo r 



since in the full Bonn model no explicit rj contribution was required by the NN data j44|. 
Furthermore, the value of the 7r- r?-m ixing matrix element is uncertain. With our particular 
choice we get^ = -3.5 x 10-^ Using g^i/An = 3.68 and {tt^\H\7]) = -5900 MeV^, as 
done in Ref. gives Pi = —1.2 x 10"^. 

One important issue is the relative sign of this contribution, which is apparently not 
determined experimentally. The sign given above is consistent with SU{3) x SU{3) chiral 
perturbation theory, which can be formulated in terms of a pseudoscalar octet tTq and a 
baryon octet. The sign of the tcs-tts niixing is, in leading order, fixed by — m^. The 
interactions of vts and vrg with the nucleon are determined by the standard weak couplings 
D and F, which are fixed in weak decays. With our definitions of g^, g>n, and jSi given above 



and the values of D and F given, e.g., in Ref. 45|, we find g^ > if we define ir^ = tt^, 
gri > and {n^\H\ri) < 0, so that Pi < 0. This conclusion holds, as it should, regardless of 
the sign definition of t], that is, whether one takes rj as vrg or — vrg. 

Fig. I2^b) represents the process where a CSB one-body operator produces a charged 
pion which then changes into a neutral pion as it re-scatters on another nucleon via the CI 
Weinberg- Tomozawa term. This contribution is small in dd avr^, since the isospin cou- 
plings force the pion exchange to occur inside one of the deuterons. This is a situation very 
similar to the seagull CSB terms, which was discussed above, but with a smaller coefficient. 
Note that a similar diagram where the exchanged pion is neutral is also small, since the on- 
shell tt'^N —>■ TC^N amplitude receives contributions only at one order higher than that from 
the Weinberg- Tomozawa term. Since the operator in Fig. I2td) is a relativistic correction to 
the leading order pion rescattering, it has exactly the same spin-isospin structure (except its 
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last term) as can be seen in Eq. Thus its first few terms are also confined to in-deuteron 
exchanges and since they are already suppressed by two orders (fcj/M)^ ~ ra-^jM, these 
terms are negligible. The last 5M /5M term has an extra Pauli spin matrix and can possi- 
bly be important since this may allow for momentum sharing. However, this term always 
includes the momentum of a final nucleon, which is very small near the pion threshold and 
this NNLO amplitude is likely to be suppressed as well. We will not consider these operators 
any further. 

The pion loops in Fig. Efe-k) represent long-range, non-analytic contributions as well as 
short-range, analytic effects. The latter cannot be separated from the short-range contribu- 
tions of Fig. 0fb), originating from a four-nucleon-pion CSB contact interaction. In this first 
study, we limit ourselves to an estimate of these effects via resonance saturation from various 
heavy- meson exchange currents HMECs — see Fig. [7|^b). In the case of the pp —>■ ppn^ reac- 
tion, heavy-meson exchanges involving the creation of a nucleon-anti-nucleon pair (z-graphs) 
were shown to be important for the total (CI) cross section near threshold These 



exchanges correspond to contact interactions in the EFT [2iJ, 21 1. Here, we include the anal- 
ogous CSB interactions where CSB occurs in the pion emission or in the meson exchange. 

The HME two-body operators are derived directly from a low-energy reduction of the 
Feynman rules for the HME Lagrangian Eq. (jSj). This gives the cr-meson-exchange two- 
body operator 

a = A^lJ:^^■iKf^J+f^J^^), (12) 

= AttM r,, ' ^^^^ 

where only the symmetrized recoil term has been used. Note that the sum is over i ^ j 
rather than i < j. 

The a;-exchange two-body operator is 

= -All ^ I*^^' ■ ^^^^^ + ^^^^^ + K^^x^,) ■ mJTj - JTM , (14) 

- if!?— ■ 

Note this has an overall minus sign and crj instead of CTj compared to Ocr. Finally there is 
a new term involving the momentum transferred to nucleon j. 
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The p-exchange two-body operator is 

fp. = . (17) 

The p HMEC is of order of the small vector times the large tensor coupling constant and 
has no contributions of order of the tensor coupling squared. 
The p-a;-mixing two-body operator is 

O,-. = -Ap-4 + -rfr^)^, ■ (Kf^ + /^k.) (18) 

+ T^rfil + C,)](^.X^,) ■ (k;/^- - i^-k,)}, (19) 
A _ _9^u^( {p\H\u;) \ 

2 

fP^ _ 9p9uj ^g; („-morii „-m^rii\ /qi \ 

AirMn.ml-ml^ ^ ' 

where the p-uj mixing is given by {p\H\uj) = —4300 MeV^ |i4Qi]. A somewhat smaller number 
{{p\H\uj) = —3500 ± 300 MeV^) was obtained in a more recent analysis 48|. The isospin- 
independent part of this p-uJ operator is only of the order of the small p vector coupling. 
Note, however, that there is a t^Tj term that involves the large p tensor coupling. 

At momenta much smaller than the heavy-meson masses these HMECs are equivalent to 
short-range pion-two-nucleon contact interactions, with specific values for the LECs. For 
example, the a mechanism [Eq. (|T!^ ] goes into the interaction shown in Eq. with 71 
given by jSig'^/ {Anm'^M) , which is consistent with the naive dimensional estimates. 

In addition, we need to consider contributions from soft photons. There is a Coulomb 
interaction and a magnetic interaction (Fig. O]), and a three-body term [Fig. Efa)]. As a first 
estimate, we shall compute the lowest order two-body diagram with a photon. This appears 
at NNLO and is shown in Fig. Et^a). 

The soft photon exchange gives a structure very similar to that of p^-uj mixing: 

+ A„ + (1 + AOrMKo-iXo-,) . (kj/ij - /Jk,)}. (22) 
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Note that the structure of this term is a consequence of gauge invariance, and this is why 
no new unknown parameters are introduced. 

III. SIMPLIFIED MODEL 

The interferences and relative importance of the CSB amphtudes of the previous section 
can be estimated in a simphfied model, using a plane-wave approximation and the simplest 
possible d and a bound-state wave functions, those of a Gaussian form. A Feynman diagram 
for this model can be drawn as in Fig. |S| Assuming spatially-symmetric bound-state wave 
functions, the invariant amplitude is given by 

M = J d^rd^pid% {A\0\DD), (25) 
\A) = ^^^^{r,puP2)\a), (26) 
\DD) = v^<l>rf(pi)<l>rf(p2)Mrf), (27) 

where and are the spatial parts of the a-particle and deuteron bound-state wave 
functions, and s = AE'j is the total cm. energy squared. The ket vectors \a) and \dd) 
contain the fully anti-symmetrized spin and isospin wave functions. Because of the symmetry 
requirements, the plane- wave dd scattering wave function is included in \dd) as given by 
Eqs. fl34|) and (jH^ below. The invariant amplitude can then be written as 

M = ^2E^s J rfWpirfV2^i(r,Pi,P2)(a|OMrf)$d(pi)$d(p2), (28) 

where {a\0\dd) contains all the spin-isospin couplings of the nucleons and the pion produc- 
tion operator O. 

The wave functions are expressed in terms of the (2+2) Jacobian coordinates 

R = -(ri + r2 + rs + r4) (= in cm.), 

r = -(ri + ra - rg - r4), 
Pi = ri - r2, 

P2 = r3-r4, (29) 
with the corresponding momenta 

K = ki + ka + ka + k4 (= in cm.), 
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FIG. 8: Feynman diagram for pion production in the dd — > an^ reaction, indicating the labeling 
of nucleons and defining basic kinematic variables. 



^(ki + ka - kg - k4) = ^(pi - ps) 
1 



p m cm. 



1^1 = -(ki - ks), 
1^2 = ^(ks - k4). 



(30) 



defined so that X]j kj ■ = K ■ R + k ■ r + ki ■ + ^2 ■ p2- The Jacobians are equal to unity 
in both representations. 

The Gaussian functions that represent the ground state wave functions are explicitly 
expressed in these coordinates using I]j<i(^j ^j)^ — + 2p^ + 2pl; 



1 



exp 
exp 



(31) 
(32) 



^3/4^3/2 2p2' ^ 

where the parameter values a = 2.77 fm and (3 = 3.189 fm are derived from measured a 
and d rms point radii; (r^)^/^ = 1.47 fm and {rf)-^^^ = 1.953 fm |49|. 

Since we have assumed that the orbital parts of the wave functions are symmetric under 
the exchange of any pair of nucleons, we may define the initial- and final-state spin-isospin 
wave functions as 

\a) = {((1, 2)i, (3, 4)Oo [[1, 2]o, [3, 4]o]o - ((1, 2)o, (3, 4)o)o [[1, 2]i, [3, 4]i] J , (33) 



\dd) = -^(l-P23-i^24)|«34), 



\di2ds,) = ((l,2)i,(3,4)i)5[[l,2]o,[3,4]o]o^(e*P-^ + (-)^e-P-^) 



(34) 
(35) 



where {i,j)s {[hjW) are the spin (isospin) Clebsch-Gordan couphngs, with magnetic quan- 
tum numbers suppressed, for nucleons, or nucleon pairs, i and j coupling to spin s (isospin 
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T). Here, Pij is the permutation operator of the indicated nucleons. The symmetry require- 
ments for the exchange of the deuterons are represented by the (orbital-angular-momentum 
dependent) combination of plane waves in Eq. (j35|) . with p as the deuteron relative momen- 
tum. Even though the expression for the a state seems to single out a (12) + (34) configu- 
ration, it is indeed fully anti-symmetric in all indices. This particular form is used because 
it closely matches the form of the initial-state wave functions, simplifying the evaluation of 
the spin-isospin summations in the matrix element. The dd wave function can in practice 
be simplified to 

\dd) = v^((l,2)i,(3,4)i)5[[l,2]o,[3,4]o]o(2L + l)z^ji(pr)^L(p-f), (36) 

since each of the three terms in Eq. ()34j) gives identical contributions to the matrix element, 
and e^^'^ + {—)^e~^^'^ reduces to 2(2L + j ^{pr) Pl{p ■ r) for any particular partial wave. 

We may obtain selection rules for the CSB amplitudes that can contribute by comparing 
this expression for the deuterons with the a-particle wave function. It is clear that matching 
the first term of |a) involves no nucleon spin or isospin fiips, but to match the second 
term, the spin and isospin of two nucleons (one from each deuteron) need to be flipped 
simultaneously. Of course, the overall spin has to change in both cases. 

In an explicit and straightforward representation, the above spin-isospin wave functions 
can be written as 



\a) 



TTii + iitr -^(Titi + itit + TiiT + itti) 



2V6 I 

(Tin + itit -TiiT -itti) 



ppnn + nnpp {pnpn + npnp + pnnp + nppn) 

2 



{pnpn + npnp — pnnp — nppn) 

(37) 



Mi2'^34lO) = ■^^--j={]]i[ - [[]]){pnpn + npnp - pnnp - nppn) ?,iji{pr)Pi^^ (38) 

where the arrows indicate spin projections and p/n proton and neutron isospin states. Note 
that for the dd state, only the spin-1, ms = state is given. These expressions can then be 
used together with the Pauli matrices of the pion production amplitudes to find the formulas 
for the matrix elements. 

In the normalization used here, the spin-averaged cross section (for s-wave pions) is given 

by 

167r.p9p^' ' ^ ^ 
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where the summation is over the deuteron polarizations. 

The CSB operators can now be evaluated in this model and studied in more detail. We 
will start with the simplest operator (i.e., the one-body term) and use its matrix element as 
a reference point for the values of the other amplitudes. 



A. One-body operator 

The one-body amplitude is strongly favored by the symmetries of initial and final states 
because all of the nucleons contribute coherently to the cross section. However, it does not 
provide momentum sharing between the deuterons and is hence dependent on the shape 
of the high-momentum tail of the ct-particle wave function. The matrix element for this 
operator is 

Ml = -i^p.(eiXe2)4>Vi, (40) 

Wi = ^2Ks j £r£p,£p2KMpr)^i^2: (41) 

where ei and 62 are the polarization vectors of the initial deuterons and the factor of 4 arises 
from the sum over all nucleons. Thus the spin- momentum structure of the ^PqS partial wave 
has been separated from the dimensionless form factor Wi. For Gaussian wave functions 
this matrix element and the corresponding cross section can be calculated analytically. They 
are 

= -z-4p.(.,xe2) ^^,^^^,^3 exp(^-^J, (42) 

where the exponential stems from the Fourier transform of the a-particle wave function and 
reflects the dependence on its high-momentum tail. We will use this one-body estimate as 
the benchmark for the calculations of more complicated amplitudes, and also for the full 
calculation using realistic wave functions. 



B. Meson-exchange operators 

Although the seagull amplitude is leading order in xPT, it is suppressed in our plane 
wave treatment of the dd — > Q;7r° reaction because of the combination of two r matrices 
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and one a matrix, which gives a poor match of the initial and final states in our simplified 
model. Thus, the pion exchange is allowed only between nucleons from the same deuteron, 
forbidding an advantageous momentum sharing between the deuterons. In addition, the 
qj term vanishes, leaving the recoil term (co'/M)kj as the only contribution. This term is 
NNLO. The pion-exchange matrix element is 

= -?^p-(eiXe2)4>V^, (44) 
= >Vi^(/r2), (45) 



(5M - \5M 

{fi2) = 1^/2^1^ / d'rd'p,d'p2Kfl2Upr)<^i<^2. (47) 



The matrix element for cr-meson exchange is 

= -z^p-(eiXe2)4W., (48) 

= / d\d'p,d'p,Kl Y: (- V. /5 + /5 V.) SMpr)p . f$i$2 (49) 

= Wi ((/r^) + (Ks) + (fis)) , (50) 

where fijjo{pr) has been replaced by ^(— Vi fij + fij Vi)ji{pr) because of the symmetriza- 
tion in Eq. (fT^ . The {ffj) are defined as the averages 

ifu) = ^^2^1^ / d'rd'p,d'p,^*J^Mpr)^i^2, 

(/fs) = :^^^^P- Jd\d'p^d'p2KfU^r + W,,-Wp,)3n{pr)p-r^i^2, 

ifu) = ^^^^P ■ / d'rd'p,d'p2K (- V. - Vp, + Vp.) /r33ji(pr)p •f<|.i$2,(51) 

where x could be any of the heavy mesons. This separates the in-deuteron exchanges /12 
from exchanges between nucleons from different deuterons (/13 and f[^). Exchanges between 
other pairs of nucleons can be reduced to these two because of the symmetries of the dd and 
a wave functions. Some of the necessary integrals are presented in the Appendix. 
The cj-exchange form factor is given by 

= wi(-(/r2)+2(/r3)), (52) 
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with the same external spin factors as in Eq. ()48p and with averages defined as in Eq. ()5ip. 
Due to cancellations, the {fi^') term cannot contribute. 

The p-exchange amplitude has an isospin factor ■ tj that makes it possible to match 
the initial-state deuterons with the second term of the a-particle wave function [Eq. (j33|) ]. 
giving sizable exchanges between the deuterons. The large ratio of the p tensor to vector 
couplings enhances this amphtude despite the small value of the p coupling constant. The 
p form factor is 

= Wi3 [(/f^) + (1 + C,)i{fQ - if?,'))] . (53) 

Similarly, the p-tu-mixing amplitude has a t^Tj term that involves the large p tensor 
coupling, which gives a large contribution to CSB for dd avr^ and possibly also for 
np (ivr^. The form factor for dd air^ is 

W,-. = Wi^[(3 + C,)(A7) - (1 + C,){f?-% (54) 

which follows immediately from the expressions for the p and uj exchanges. The /12 term 
vanishes since the 1 + r| term of Eq. (PT|) gives zero when acting on a deuteron and the 
(Ti X (T2 term vanishes due to the spin-couplings in the wave functions. 
The photon exchange contribution is 

= 1^1x7'^^ ^ ^^"^^^"'^ " ^ 2A,)(AV)], (55) 

where again the fu term vanishes. In the simplified model the photon exchange only oc- 
curs between pairs of protons, thus not benefiting from the coherence the other amplitudes 
experience. However, the relatively large coupling makes this amplitude important. 

Thus all HMECs have contributions to dd an^. While there are some cancellations of 
the in-deuteron (/12) and the derivative (^3) exchange terms between the different heavy 
mesons, all the /13 contributions are of the same sign. Since in all cases (/13) is much larger 
than (/12) and (/(3) (see Table HTl below), they dominate the matrix element and the cross 
section, adding coherently with each other and also with the one-body and pion-exchange 
terms. Also the photon graph is of the same sign. Thus the internal spin-isospin symmetries 
of the dd : a system used here strongly favor the one-body and meson-exchange amplitudes if 
the plane wave approximation is used. This will be demonstrated quantitatively in Sec. IIVI 
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TABLE II: The Yukawa averages (see text for definitions) evaluated for dd — > an^ at the two lUCF 
energies, using Gaussian wave functions. 







Td = 228.5 MeV 






= 231 8 MeV 






\J 1^/ 


\J u/ 


\J 13/ 


\J 1^/ 


\J Li/ 


\J 13/ 


TT 


0.0172 








0.0172 








a 


0.0292 


0.470 


0.0220 


0.0292 


0.490 


0.0229 




0.0228 


0.395 


0.0106 


0.0228 


0.412 


0.0111 


P 


0.00095 


0.0165 


0.00046 


0.00095 


0.0172 


0.00048 


p-uj 


0.00445 


0.0704 


0.00377 


0.00445 


0.0734 


0.00393 


7 


0.00073 


0.0032 


0.00109 


0.00073 


0.0033 


0.00111 



IV. RESULTS 



The matrix elements of the previous section are evaluated numerically using the simplified 
(Gaussian) wave functions. At most a double integration with a separate single integral was 
needed, which was carried out using standard Gauss-Legendre techniques. Explicit formulas 
for the integrals are presented in the Appendix. 

The Yukawa averages are tabulated in Table |n] for each operator and both energies 
relevant to the lUCF experiment. The (/12) contributions are the same at both energies 
since the in-deuteron exchange is independent of the energy in our plane-wave model. For 
completeness the (/12) values for p-uj and photon exchanges are given, even though these, 
as discussed above, do not contribute to the matrix element in our simplified model. Note 
that even though the p-exchange, p-cu-mixing, and photon integrals are much smaller than 
those for the other meson exchanges, they will be multiplied by large constant factors in the 
definitions of the matrix elements, Eqs. ()53 p - ()55|) . This drastically increases their relative 
importance. 

The matrix elements and cross sections calculated from these averages are given in Ta- 
ble mil individually for each amplitude and as a grand total. For comparison purposes, 
the experimental cross sections are included and the matrix elements are given relative to 
the one-body matrix element. All the heavy meson exchanges are of the same order as the 
one-body term, with the p-u mixing being the largest. Adding all amplitudes gives a total 
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TABLE III: Matrix elements and cross sections evaluated for dd an at the two lUCF energies. 
The matrix elements are given relative to the one-body matrix element, with the relevant CSB 
mechanism indicated. The experimental cross sections are also given. 





M(228 5) 


M(231 8) 


crf228 5) 


C7f231 8) 






\M^] 
[•' ' IJ 


fob! 


lir"~'\ 


vr (SM SM) 


0.128 


0.128 


0.011 


0.014 


1 ( TT-tl 1 


1 

X 




688 


n 869 


a [tt-t]) 


0.522 


0.543 


0.187 


0.256 


UJ {ir-rj) 


0.766 


0.801 


0.404 


0.557 


P (t^-v) 


0.344 


0.359 


0.082 


0.112 


p-u (p-uj) 


1.546 


1.612 


1.645 


2.256 


7 (el. -mag.) 


1.469 


1.517 


1.486 


1.999 


total 


5.78 


5.96 


23.0 


30.8 


Exp. \U} 






12.7 ±2.2 


15.1 ±3.1 



matrix element that is almost six times that of the one-body, increasing the cross section 
from a meager 0.69 pb at Ta = 228.5 MeV (0.87 pb at 231.8 MeV) to 23 pb (31 pb). This is 
of the same order as the lUCF data cr = 12.7 ±2.2 pb and 15.1 ±3.1 pb. Note that the cross 
section is not strictly linear in the pion momentum — the momentum transfer in the wave 
functions introduces a dependence on the deuteron momentum, which modifies the linearity, 
at least in this simplified model. For example, Eq. ()43|) for the one body term contains the 
square of the deuteron momentum in an exponential. 

The relative proportions of the pion-exchange {SM — ^6M), photon exchange, p-uj-, 
and TT-r^-mixing (sum of one-body and HMEC) contributions to the matrix element are 
roughly 7r:7:p-co':7r-?7=l:ll:12:21. Thus the formally leading seagull terms make up only 
about 2% of the total matrix element. The total cross section can be expressed in terms of 
the relative contributions of the different CSB mechanisms such that the dependences on 
the corresponding parameters are made explicit: 

(t(228.5 MeV) = (23.0 pb) (0.254 ± 0.0188—^^— ± 0.0034- 



2.03 MeV -0.74 MeV 
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+ 156 ^''^^ ^'^'^'''"^ + 26S ^'^^ - ^^'^1^°^ ^ f56) 
V47r • 0.51 (-4200 MeV^) ' 47rV0.43 ■ 10.6 (-4300 MeV2) j " ^ ' 

Here the numerical coefficients are the fractions of the matrix element belonging to each of 
the considered mechanisms, assuming our choice of parameter values. The various terms are 
normalized, as indicated, to these values. The photon exchange diagram is represented by 
just a number, since its parameters are well-known. Note that the second and third terms 
are constrained by the neutron-proton mass difference [Eq. (jS))]. The ?7-7r°-mixing term can 
be further separated to show the relative contribution of the various HMEs. Thus, 

^ = (2.18 v^) (o.380 + 0.198j^ + 0.291j^ + 0.131^) . (57) 

where a-j^n is the cross section from ti-t] contributions alone and the ffist number in the 
second parenthesis is the one-body contribution. At the higher lUCF energy the relative 
weights of the different contributions in Eqs. (|56p and (j57|) remain more or less the same, 
with only minor changes. The sensitivity of the cross section calculation to a different choice 
of couplings can easily be found from these two formulas. For example, using instead the 
large g^/4:TT = 3.68 and {ri\H\7c^) = —5900 MeV^ as in Ref. 0|, the cross section increases 
from 23 to 118 pb (31 to 158 pb at 231.8 MeV). 



V. DISCUSSION 



Our simplified model keeps a complete treatment of the dominant pieces of the spin- 
isospin couplings in the bound state wave functions, even though it ignores some dynamics 
of the dd —>■ avr" reaction and the distortion of the initial state. As a result, the symmetries 
of the bound state wave functions allow us to determine the CSB amplitudes that are 
guaranteed to be important for a full calculation. 

This treatment shows that the LO 7r-rescattering term (from the the seagull interactions) 
is suppressed because of a poor overlap with the initial- and final-state wave functions. 
Photon loops, at NLO, vanish due to symmetries and cancellations, while a three-body con- 
tribution might survive, but has not yet been calculated. On the other hand, the NNLO 
one-body amplitude and N^LO heavy-meson exchanges are strongly favored, adding coher- 
ently with each other. Their dominance would be even more spectacular if a larger value 
for the rjNN coupling were used. Also the p-cu-mixing and photon-exchange terms are 
important and enter at the same level as the one-body and HMEC terms. 
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We note that our analysis assumes that pions are produced in s-waves, as one would 
expect for a near-threshold reaction. This is supported b y th e lUCF experiment, where the 
energy dependence is consistent with s-wave production jl7|. 

If our simple wave functions and the plane wave approximation are used, then, within 
the resonance saturation picture, the dominant CSB mechanisms for the dd — > avr^ reaction 
are identified with tt-t] mixing (one-body enhanced by HMECs), followed by p-uo mixing 
and photon exchange, and finally a small contribution from pion rescattering (related to the 
neutron-proton mass difference). 

The coherent sum of these contributions leads to a cross section of the same order of 
magnitude as the observed one 17[, but more needs to be done before making a detailed 
assessment of the quality of the agreement between theory and experiment. 

It is likely that the relative importance of these amplitudes will be shifted once realistic 
wave functions are used. Preliminary calculations suggest that the one-body term can in 
fact be enhanced by as much as a factor of three or four with a realistic a-particle wave 
function. This is expected, since this amplitude is sensitive to the high-momentum tail of 
the wave function, which is very small for a Gaussian. Preliminary estimates also show that 
spin-dependent initial state interactions enhance the pion rescattering contribution. We 
stress that the HMECs are less sensitive to the a-particle wave function and should remain 
crucial for the interpretation of CSB in the dd —>■ avr" reaction. 

A full model calculation, using realistic bound-state and drf-scattering wave functions is 
needed in order to have a clear understanding of the CSB mechanisms behind the dd avr^ 
reaction. Furthermore, the effects of CSB in the wave functions and some diagrams ignored 
here, such as the long-range part of the various NNLO pion loops and the N^LO recoil part 
of the A-excitation term, should be included. In particular, it is necessary to include the 
non-z-graph part of photon exchange, e.g., the Coulomb interaction in the initial and final 
states (Fig. E)). Such an investigation is currently in progress and will be reported later. 
The general conclusions and insights from the present paper provide important guidelines 
for that work. 

We note a very interesting parallel between the dd an^ process considered here, and 
the reaction pp ppn . In both formally-leading diagram is suppressed and the 

sub-leading diagrams are crucial to explain the cross section. Despite several serious efforts 
that have yielded substantial insights into the various A^A^ NNtt systems, the pp — ppn^ 
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reaction is still not completely understood, especially regarding spin observables [s^ . 

Higher-order interactions — such as the heavy-meson-exchange terms, which could in- 
crease the role of TT-rj mixing — might help improve the agreement between the TRIUMF 
result for A{^(np dir^) [l6| and theoretical estimates based on reasonable values for SM 
and 6M It is thus necessary that a future calculation of Af^(np dn'^) includes these 
higher-order terms. The three reactions pp — » ppn^, np dn^, and dd an^ provide 
important testing grounds for any pion-production model that intends to include effects 
beyond- leading-order. 
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APPENDIX: EXPLICIT EXPRESSIONS FOR YUKAWA AVERAGES 

The averages of the different Yukawa factors of Eq. can be reduced to at most two- 
dimensional integrals, using the Gaussian wave functions of Sec. IIIII The angular and one 
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of the radial integrals can be carried out analytically, resulting in the explicit formulas 
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j dpp^e p j drr^jo{pr)e 



(A.l) 



where ri3 = |ri — r3| and 1/7^ = + 1/(2/3^). 



[1] G. A. Miller, B. M. K. Nefkens, and I. Slaus, Phys. Rep. 194, 1 (1990). 

[2] S. Weinberg, Trans. N.Y. Acad. Sci. 38, 185 (1977). 

[3] H. Lcutwylcr, Phys. Lett. B378, 313 (1996). 

[4] L. M. Barkov et al, Nucl. Phys. B 256, 365 (1985). 

[5] J. A. Nolen and J. P. Schiffer, Ann. Rev. Nucl. Part. Sci. 19, 471 (1969). 

[6] B. Gabioud et al, Phys. Rev. Lett 42, 1508 (1979); Nucl. Phys. A420, 496 (1984). 

[7] R. Abegg et al, Phys. Rev. Lett. 56, 2571 (1986); R. Abegg et al, Phys. Rev. D 39, 2464 
(1989); S. E. Vigdor et al, Phys. Rev. C 46, 410 (1992); R. Abegg et al, Phys. Rev. C 57, 
2126 (1998); G. A. Miller, A. W. Thomas, and A. G. Wilhams, Phys. Rev. Lett. 56, 2567 
(1986); A. G. Williams, A. W. Thomas, and G. A. Miller, Phys. Rev. C 36, 1956 (1987). 

[8] N. Fettes and U.-G. Meifiner, Nucl. Phys. A 693, 693 (2001); and references therein. 

[9] L. I. Lapidus, Zh. Eksp. Tcor. Fiz. 31, 865 (1956) [Sov. Phys. JETP 4, 740 (1957)]. 
[10] N. E. Booth, O. Chamberlain, and E. H. Rogers, Buh. Am. Phys. Soc. 4, 446 (1959); Nuovo 
Cimento 19, 853 (1961); Y. K. Akimov, O. V. Savchenko, and L. M. Soroko, Zh. Eksp. Teor. 
Fiz. 41, 708 (1961) [Sov. Phys. JETP 14, 512 (1962)]; J. A. Pohier and M. Pripstein, Phys. 
Rev. 122, 1917 (1961); 130, 1171 (1963). 

33 



[11] J. Banaigs et al, Phys. Lett. B53, 390 (1974); Phys. Rev. Lett. 58, 1922 (1987). 
[12] L. Goldzahl. J. Banaigs, J. Berger, F. L. Fabbri, J. Hiifner, and L. Satta, Nucl. Phys. A533, 
675 (1991). 

[13] P. Fleury and F. Plouin, "Status Report LNS Expt. 105", (1987) (unpublished). 

[14] D. Dobrokhotov, G. Faldt, A. Gardestig, and C. Wilkin, Phys. Rev. Lett. 83, 5246 (1999). 

[15] A. Gardestig, G. Faldt, and C. Wilkin, Phys. Rev. C 59, 2608 (1999). 

[16] A. K. Opper et al., Phys. Rev. Lett. 91, 212302 (2003). 

[17] E. J. Stephenson et al, Phys. Rev. Lett. 91, 142302 (2003). 

[18] V. Bernard, N. Kaiser, and U.-G. Meifiner, Int. J. Mod. Phys. E 4, 193 (1995). 

[19] P. F. Bedaque and U. van Kolck, Ann. Rev. Nucl. Part. Sci. 52, 339 (2002). 

[20] T. D. Cohen, J. L. Friar, G. A. Miller, and U. van Kolck, Phys. Rev. C 53, 2661 (1996). 

[21] U. van Kolck, G. A. Miller, and D. O. Riska, Phys. Lett. B388, 679 (1996). 

[22] C. A. da Rocha, G. A. Miller, and U. van Kolck, Phys. Rev. C 61, 034613 (2000). 

[23] C. Hanhart, U. van Kolck, and G. A. Miller, Phys. Rev. Lett. 85, 2905 (2000). 

[24] C. Hanhart and N. Kaiser, Phys. Rev. C 66, 054005 (2002). 

[25] B. Y. Park, F. Myhrer, J. R. Morones, T. Meifiner, and K. Kubodera, Phys. Rev. C 53, 1519 
(1996); T. Sato, T.-S. H. Lee, F. Myhrer, and K. Kubodera, Phys. Rev. C 56, 1246 (1997); 
C. Hanhart, J. Haidcnbaucr, M. Hoffmann, U.-G. Meifiner, and J. Spcth, Phys. Lett. B424, 8 
(1998); V. Dmitrasinovic, K. Kubodera, F. Myhrer, and T. Sato, Phys. Lett. B465, 43 (1999); 
E. Gedalin, A. Moalem, and L. Razdolskaya, Phys. Rev. C 60, 031001 (1999); V. Bernard, 
N. Kaiser, and U.-G. Meifiner, Eur. Phys. J. A4, 259 (1999); S. Ando, T.-S. Park, and D.- 
P. Min, Phys. Lett. B509, 253 (2001). 

[26] U. van Kolck, Few-Body Syst. Suppl. 9, 444 (1995); Ph.D. Dissertation, U. of Texas (1993). 

[27] S. Weinberg, in Chiral Dynamics: Theory and Experiment, edited by A. M. Bernstein and 
B. R. Holstein (Springer- Verlag, 1995). 

[28] G. Ecker, J. Gasser, A. Rich, and E. de Rafael, Nucl. Phys. B 321, 311 (1989). 

[29] E. Epelbaum, U. G. Meifiner, W. Glockle, and C. Elster, Phys. Rev. C 65, 044001 (2002). 

[30] U. van Kolck, J. L. Friar, and T. Goldman, Phys. Lett. B371, 169 (1996). 

[31] C. Y. Cheung, E. M. Henley, and G. A. Miller, Nucl. Phys. A 305, 342 (1978). 

[32] J. A. Niskanen, Few-Body Syst. 26, 241 (1999). 

[33] U. van Kolck, J. A. Niskanen, and G. A. Miller, Phys. Lett. B493, 65 (2000). 

34 



[34] S. A. Coon and B. M. Freedom, Phys. Rev. C 33, 605 (1986); see also a calculation applicable 

to the energy region of the (3,3) resonance: C. Y. Cheung, Phys. Lett. B119, 47 (1982). 
[35] J. Gasser and H. Leutwyler, Phys. Rep. 87, 77 (1982). 

[36] V. G. J. Stoks, R. A. M. Klomp, M. C. M. Rentmeester, and J. J. de Swart, Phys. Rev. C 48, 

792 (1993); U. van Kolck, M. C. M. Rentmeester, J. L. Friar, T. Goldman, and J. J. de Swart, 

Phys. Rev. Lett. 80, 4386 (1998). 
[37] E. Epelbaum, A. Nogga, W. Glockle, H. Kamada, U.-G. Meifiner, and H. Witala, Phys. Rev. 

C 66, 064001 (2002); Eur. Phys. J. A15, 543 (2002). 
[38] All OBE parameters are adapted from T. E. O. Ericson and W. Weise, Pions and Nuclei, 

(Oxford University Press, New York, 1988). 
[39] C. Hanhart, G. A. Miller, F. Myhrer, T. Sato, and U. van Kolck, Phys. Rev. C 63, 044002 

(2001). 

[40] S. A. Coon and M. D. Scadron, Phys. Rev. C 51, 2923 (1995). 
[41] L. Tiator et al, Nucl. Phys. A580, 455 (1994). 
[42] O. Dumbrajs et al, Nucl. Phys. B216, 277 (1983). 

[43] R. Machleidt, C. Elster, and K. Holinde, Phys. Rep. 149, 1 (1987); R. Machleidt, Adv. Nucl. 

Phys. 19, 189 (1989). 
[44] R. Machleidt, Phys. Rev. C 63, 024001 (2001). 
[45] R. P. Springer, Phys. Lett. B461, 167 (1999). 
[46] T.-S. H. Lee and D. O. Riska, Phys. Rev. Lett. 70, 2237 (1993). 
[47] C. J. Horowitz, H. O. Meyer, and D. K. Griegel, Phys. Rev. C 49, 1337 (1994). 
[48] S. Gardner and H. B. O'Connell, Phys. Rev. D 57, 2716 (1998); 62, 019903(E) (2000). 
[49] J. Carlson and R. Schiavilla, Rev. Mod. Phys. 70, 743 (1998). 
[50] H. O. Meyer et al, Phys. Rev. C 63, 064002 (2001); C. Hanhart, hep-ph/0311341. 
[51] This assignment has no deeper reason than the numerical proximity of \/m^M = 2.6m^ and 

Ma - M = 2.1m^. 

[52] Apart from the pion mass itself, the effects of explicit chiral symmetry breaking are typically 
sub-leading. As a consequence, quark-mass isospin breaking is in g eneral smaller than it would 
be expected on the basis of the quark mass splitting, e ~ 1/3 |26| |. 

[53] Note that the leading ^ttNN vertex — the so-called Kroll-Rudermann term — couples to the 
pion charge and therefore does not directly contribute to neutral-pion production. 



35 



